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Abstract 

We describe the quantum state of a Bose-Einstein condensate at zero tempera- 
ture. By evaluating the Q-function we show that the ground state of Bose-Einstein 
condensate under the Hartree approximation is squeezed. We find that multimode 
Schrodinger cat states are generated as the condensate evolves in a ballistic expan- 
sion. 
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1 Introduction 



The realisation of Bose-Einstein condensation in dilute gases has made it pos- 
sible to study the dynamics of quantum fields directly [1,2]. In this paper, we 
address one of the most intriguing issues surrounding a Bose-Einstein con- 
densate (BEC), namely the nature of its quantum state. Number states might 
seem to be the natural choice given that atoms cannot be created or destroyed 
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at these energies. However, it is well known that, in an open environment, par- 
ticles may be added or removed from the condensate, and this implies that 
the condensate is not in a pure Fock state. On the other hand, it has been 
suggested that there should be a highly entangled state of the condensate plus 
the environment [3]. The coherence properties of the condensate have been 
demonstrated experimentally, indicating that a coherent state might be the 
most robust representation [4,5,6]. This point turns out to be very important, 
as it is well-known that a coherent state propagating through an amplitude- 
dispersive medium evolves into a superposition of two coherent states [7]. The 
generation of these Schrodinger cat-hke states is considered to be a general 
property of nonlinear systems where dissipation effects are not very large. 
Recently, some schemes to create Schrodinger cat states in BECs have been 
proposed [8] and the effects of loss have also been considered. 

It has also been pointed out that the state of the condensate might be squeezed, 
given that binary collisions between atoms play a very important role in the 
description of BECs [9] and are known to give rise to squeezing [10,11,12]. We 
note that the question of describing the quantum state of a trapped BEC was 
addressed in reference [11] where a single mode approximation was considered, 
based on a symmetry-breaking picture. The authors showed that for repulsive 
interactions the state is number squeezed. In this paper, we use a multimode 
description in the Hartree approximation. In particular, we show how the 
state evolves under ballistic expansion, which is important as this is the basis 
of many experiments involving the condensates [4,5,6]. 

To study the evohition of the quantum state of a BEC, we have used the 
Heisenberg equation of motion obtained from the many-body Hamiltonian 
describing the condensate, and used the Hartree approximation, under which 
the n-atom wavefunction is written as a product of n single atom wavefunc- 
tions. This results in a nonlinear Schrodinger equation for the condensate 
wavefunction, or the well-known Gross-Pitaevskii equation (GPE). The GPE 
has proved to be a very useful tool for studying the condensate dynamics, 
and its predictions are in very good agreement with experimental results at 
low temperatures [13,14]. It is important to note that the GPE can be under- 
stood in terms of either coherent states or number states. Here, we make use 
of a number state description in order to obtain an analogue of the scenario 
presented in [15] for the case of light. We shall take "snapshots" of the wave- 
function at various times and analyse them using the Q-function, which is a 
quantum-mechanical phase space function. We also want to mention that we 
are treating the condensate in a zero temperature regime; this means that we 
are basically referring to the predominant mechanisms occurring in the ground 
state of the condensate. At finite temperature, the mechanism is affected by 
the presence of excited states and therefore the dynamics are changed [16]. 

The paper is organised as follows: in Section 2, we present the Hamiltonian 
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that describes the BEC and solve the equation of motion in the Schrodinger 
picture using the Hartree approximation. In Section 3. we calculate the Q- 
function for the condensate. The evolution of the state of the condensate in a 
baUistic expansion is then presented, and the generation of Schrodinger cat- 
like states is observed. This is relevant in studies of output couplers or atom 
lasers. Finally we conclude in Section 4. 



2 Q-Function in the Hartree approximation 

The phase space representation is useful for visualising the evolution of a 
quantum state, as the statistics of the state can easily be described using 
the quasiprobability distribution functions [17]. In this case, we have made 
use of the Q-function which has the nice property of being always positive. 
This function can be considered to describe probability densities, and has the 
further advantage of being readily measurable by quantum tomographic tech- 
niques [12,18]. The first experimental measurement of quantum states was 
reported by Smithey et. al. in 1993 [19]. Since then, several proposals to re- 
construct the quasiprobability distributions of a quantum system, such as the 
Wigner function, have been made. For instance, Kurtsiefer et. al. [20] reported 
a technique to measure the Wigner function of an ensemble of helium atoms 
in a double-slit experiment. Moreover, a proposal has been made to measure 
the quantum state of a BEC in an atomic interferometer based on Raman 
transitions [12], where the Q- function can be measured directly. The way of 
obtaining the quasiprobability function is analogous to doing unbalanced ho- 
modyne detection of an optical signal. 

The Q- function of a pure quantum state |^) is defined as 



where the two real numbers a,, and ttj describe the coherent state \a). There- 
fore, the quasiprobability distribution depends directly on these two parame- 
ters, which span the phase space. An alternative form of equation (1), leads 
the Q-function to be defined by the diagonal matrix elements of a density 
operator p in a pure coherent state \a): 

Q{ar,a.i) = -{a\p\a). (2) 



Clearly, the Q-function is characterised by being always positive and nor- 
malised to unity. Thus, it can be regarded as describing probability densities. 





TT 
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and the numbers and ctj play the role of conjugate variables, eg number 
and phase. 

We shall start by describing the system in terms of the many-body interacting 
Hamiltonian 



H = J d3r0^(r,T) 



4>{r,T) 



-i J dhdh'(i>\r, r)0^(r', r)V(r - r')0(r, r)0(r', r), 



(3) 

(4) 



where 0^(r, r) and 0(r,r) are the boson field creation and annihilation oper- 
ators respectively, and V{r — r') is the two-body interatomic potential. For a 
contact interaction potential, the normalised Heisenberg equation of motion 
for the operator (p is given by 



90 



-V^ + Vtrap{r) + 



nU< 



rap 



4>- 



(5) 



In the Schrodinger picture, equation (5) takes the form i/i^|$) = H\^). The 
solution can then be found by expanding the state vector |$) in a Fock space 



1$) = Y,A,,\n]t). 



(6) 



The solutions of the equation are then number states |n;t). In the Hartree 
approximation, it is possible to write the number states as 



\n\t) 



dx'^n{x,t)^\x)) |0), 



(7) 



where ^„ is a single particle wavefunction and 0^ corresponds to the creation 
operator of the field. 

This implies that each particle in the system experiences the same potential. 
Therefore, the single-particle function ^„ satisfies the GPE with a scaled 
nonlinearity proportional to the number of particles [15]. 

If we superimpose these states with a Poissonian distribution, the coefficient 
An takes the form 



an 



exp 



"0 



(8) 
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so that the quantum state in the Hartree approximation is expressed as 

|2 



l*) = E^exp 



n 



(9) 



The Q-function given by equation (1) can be written in terms of a reference 
coherent state \a,{^{x,t)}) as 



(10) 



We can construct the reference coherent state |q;, t)}) using the n- 

particle eigenstate for the field with envelope 'i/n{x,t), where n = |q;oP has 
the meaning of the average particle number 



= ^ (/ dx*^(x,i)0t(^)y |o). 



(11) 



Assuming Poissonian statistics, we obtain the many-body coherent state 



a,{'iffi{x,t)}) = exp ( ) exp (a / dx'ifn{x,t)4>\x)]\0). 



[12) 



Using equations (9), (10) and (12), the Q-function can be written 



Q(ar, ai, t) = exp (-|q;o|^ - |q;|^) 



X 



n=0 



ni 



J dx%{x,t)-^n{x,t)j 



.(13) 



3 Quantum state of a Bose-Einstein Condensate at zero tempera- 
ture 



Let us consider Gaussian wavefunctions to describe the ground state of the 
trapped condensate 



exp 



X 



n 



(14) 



and assume that n is large so that the difference between n — 1 and n can be 
neglected. The Q-f unction for such a case is 



Qaioir, ai) = exp (-|q;o|^ - |q;|^) 



n=2 



ni 



_2 2 \ in- 



.(15) 
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For an average number of particles n — 50 the Q-function (15) is shown in 
Figure 1. 
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Fig. 1. Q-function for a Bose condensed gas described by Gaussian functions. 



A more reahstic approximation to the ground state is given by the Thomas- 
Fermi solution 

*n = , ^ il^n - Vtrapf^ , (16) 

^Jn{n - l)Uo 



where the chemical potential is determined by the normalisation of 
The Q-function is then expressed as 



n=2 



(A*n - A«n)F ( J — , 
2 




iHn + IJ'n)^ [J—, J — 



(17) 



where F((/9, A;) and E((/9, A;) are incomplete elliptic integrals of the first and 
second kind respectively [21]. For a harmonic trap potential V\ = ^muf^^pr'^, 

the chemical potential is given by /i„ = ^Miizili^^ ^' ^]^q Q-function for the 
same average number of particles as in the previous case is plotted in Figure 
2. 



From Figures 1 and 2, it is clear that the ground state of the condensate 
is very close to a squeezed state, since the Q-function has different widths 
in the and aj-directions. On the one hand, the squeezing is generated by 
the nonlinearity inherent to the system. This can be explained in terms of a 
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Fig. 2. Q-function for a Bose condensed gas in the Thomas-Fermi approximation. 

superposition of different numbers of particles with different phases, creating 
a deviation from the classical phase. On the other hand, we know that the 
coherent state interpretation arises from the robustness of these states in the 
presence of interactions between the condensate and its environment [3]. It 
is well known that coherent states are a particular case of a more general 
class of minimum-uncertainty states, namely the squeezed states. It has been 
shown that the propagation of a coherent state in a nonlinear medium evolves 
into a superposition of two coherent states [7], while the quantum- mechanical 
superposition of two coherent states has been shown to lead to a transition 
from a Poisson distribution to a sub-Poissonian one and squeezing [22]. In 
this sense, the robust description of the quantum state still applies, bearing 
in mind that the existence of the interactions between atoms yields squeezing 
of the quantum field and, as a result, entanglement between the atoms is also 
produced [9]. We note that while a single coherent state does not evidence 
squeezing, a considerable amount is exhibited in a superposition of even just 
two coherent states. 

This mechanism of achieving squeezing is the predominant one in the case 
of the ground state, without the presence of highly populated excited states. 
When we consider the presence of a thermal cloud, we can treat the system 
using the Bogoliubov procedure and analyse the behaviour of the collective 
excitations. It has been shown that squeezed quasiparticle excitations can be 
produced when considering the generation of correlated pairs by the interac- 
tion of an excited state with the ground state [16]. Thus, we can talk about 
two different aspects of nonclassicality in the system, the predominant one for 
the ground state at zero temperature, and the prevalent one in the presence 
of excited states. 

It is of particular interest to see how the Q-function behaves once the con- 
densate is released from the trap as when the atomic cloud is allowed to 
expand ballistically [1], or in some output coupling schemes for atom lasers 



[23]. The Thomas- Fermi solution was used as an initial condition to solve the 
GPE after the trapping potential is switched off and the Q-function (f 7) was 
then calculated from the dynamical evolution of the wavef unctions. We expect 
that once the trapping potential is switched off, the condensate will rether- 
malise. Comparing the quantum state of the condensate, for a trap frequency 
Wtrap = '2.71 X 100 Hz, at 2 X 10^'^ s (shown in Figure 3) to the Q-function for a 
thermal state (Figure 4) makes it clear that rethermalisation is indeed taking 
place. 
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Fig. 3. Q-function for the dynamical evolution of the system at 2 x 10~^ s. 
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Fig. 4. Q-function for a typical thermal state 

A very important aspect that has to be considered is the finite lifetime of 
the condensate, which is at most 10-20 s. The reason behind this, is that the 
thermal cloud is not the only source of decoherence. The condensate loses 
atoms because of Rayleigh scattering, external heating, and three-body de- 
cay. Information about the quantum state of the condensate is carried by the 
atoms that escape from it so that, due to interactions of the system with the 
environment, a quantum superposition is turned into a statistical mixture. 
From the processes mentioned above, three-body recombination is the most 
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important one [24,25]. In an experiment carried out by Stamper-Kurn et. al. 
[26], the measured loss rate per atom turned out to be 4/s for = 10^ atoms. 
This rate scales as N^, where N is the number of atoms in the condensate. 
Considering a condensate with N — 10^, just one atom is lost per second; the 
decoherence time is 1 s. Bearing this fact in mind, it is perfectly sensible to ask 
whether, after opening the trap, it is still possible to see a Schrodinger cat-like 
state. When the states evolve further, a mult i- component structure develops 
as shown in Figures 5 and 6, at 5 x 10~^ s and 1 x 10~^ s, respectively. In both 
figures, the quasiprobability function demonstrates that the quantum states 
obtained correspond to Schrodinger cat states hke the ones generated in Kerr 
media [27,28]. Moreover, we also see that the cat states have a different num- 
ber of components at different times in the evolution, in a close analogy to 
the mechanism discussed by Yao for the case of a radiation field propagating 
in a nonlinear medium [29]. It has been shown that the quantum-mechanical 
superposition of states can yield a state that exhibits squeezed fluctuations 
[30,22]. This suggests that, even in the presence of losses due to three-body 
recombination, it would still be possible to observe Schrodinger cat-like states. 
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Fig. 5. Q-function for the dynamical evolution of the system at 5 x 10 ^ s. 



4 Discussion 

In summary, we have described the ground state of a Bose condensed gas 
in the Hartree approximation in a phase-space representation using the Q- 
function. We have confirmed that the ground state of the condensate is in 
a squeezed state due to the presence of interactions between the atoms. We 
have also studied the dynamical evolution of the condensate in a ballistic 
expansion. The simulations show that immediately after the condensate is 
released from the trap the system rethermalises, with the quantum state being 
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Fig. 6. Q-function for the dynamical evolution of the system at 1 x 10 ^ s. 

still close to a squeezed state. Continuing the evolution without the trapping 
potential, we obtain a multimode Schrodinger cat-like structure. In general, 
the "longevity" of these states is very small due to decoherence. However, 
the time scales of the calculations suggest that it might still be possible to 
observe them. The generation of nonclassical states by the time evolution of 
an initial coherent state under the influence of some nonlinear dynamics has 
been studied in the past for the case of radiation fields [7,27,28,29]. In the 
case of a BEC the nonlinearity arises from binary correlations between the 
atoms, producing the effective interactions between them. We point out that 
this process basically refers to the ground state, as the presence of excited 
states changes the predominant mechanisms of achieving squeezing. 

An issue that remains to be addressed is the one related to the experimental 
reconstruction of the quantum state of the condensate. We have mentioned 
that it is possible to use quantum tomographic techniques to measure the 
Q-function [12,18]. As an alternative to these techniques, a scheme using pro- 
jection synthesis was proposed by Baseia and co-workers [31]. The scheme 
follows the idea of measuring the probabihty distribution of a state |/) to 
have an observable 9 



where the expectation value is replaced with the probability density by setting 
Cn to zero for n > N, with large N. The probability density is proportional 

to the expectation value of the projector tt = k,\6){6\ where k is a positive 
constant. If the field is in a mixed state p, equation (18) can be written as 



1 

{f\9){9\f)^P^{0)^- ^CneM-inO) 



(18) 



n=0 




(19) 
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It has been shown [32] that by specifying a suitable reference state \B)b in an 
experimental setup where a state \ilj)a is coherently mixed in a beam splitter 
with the reference state \B)h, we can obtain the required probability (19). 

In the case of the quasiprobability function Q, we would require a projector 
TT — K\a) {a\ so that its substitution in equation (19) yields 

Pp{a) = 'IV(p7r) = {a\p\a) = Q{ar, a^) (20) 



The form of the reference state that leads to the required projector is given 
by [31] 

\B),^Cj2h\k)b, 

k=0 



where C is a normalisation constant and 

-N ( V^'"^ 



It is important to mention that this scheme involves projections onto truncated 
coherent states, which is a good approximation for large N . In other words, it 
works well for the case N ':$> n, where n is the average particle number of the 
synthesised coherent state. The construction of the reference state (21) would 
be the main experimental challenge in this scheme. However, one might able 
to generate such a state using the scheme proposed in [33] where the authors 
show that it is possible to prepare an arbitrary (finite) superposition of Fock 
states by applying a well-defined succession of the displacement and creation 
operators. Experimentally a sequence of beam splitters at which conditional 
measurements are performed is required. In the case of BEC, this may be 
achieved by using repeated interference of condensates with different hyperfine 
states, similar to the Gedanken experiment proposed by Castin and Dalibard 
in reference [34] where atoms leak from two trapped condensates and are 
detected in the output channels of a 50-50 beam splitter. Exact procedure for 
such state preparation is the subject of our future work 

This work has been supported by CONACyT, EPSRC, the Royal Commission 
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